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We use a sample of 384 million BB events collected with the BABAR detector at the PEP-II e~^e~ 
collider to study angular distributions in the rare decays B -^ K*£'^£~ , where i^£^ is either e'^e^ or 
fi'^ fj,~ . For low dilepton invariant masses, mu < 2.5 GeV/c , we measure a lepton forward-backward 
asymmetry Afb = 0.24lo;23 ± 0.05 and K* longitudinal polarization Fl = 0.35 ± 0.16 ± 0.04. For 
mu > 3.2 GeV/c^ we measure Afb = 0.76lo:32 ± 0.07 and Fl = 0.71 

PACS numbers: 13.20.He 



The decays B -^ K*£'^l~, where K* ^ Kn and 
£~^£~ is cither an e~^e~ or (J.^ fJ-~ pair, arise from flavor- 
changing neutral currents (FCNC), which arc forbidden 
at tree level in the Standard Model (SM). The lowest- 
order SM processes contributing to these decays are the 
photon or Z penguin and the W^W~ box diagrams 
shown in Fig. [1] The amplitudes can be expressed in 
terms of effective Wilson coefficients for the electromag- 
netic penguin, Cj*^*^, and the vector and axial-vector elec- 
troweak contributions, Cg" and CJq respectively, arising 
from the interference of the Z penguin and W^W~ box 
diagrams [l|. The angular distributions in these decays 
as a function of dilepton mass squared q^ 



"e+e- 



are 



sensitive to many possible new physics contributions [2]- 
We describe measurements of the distribution of the 
angle 6k between the K and the B directions in the K* 
rest frame. A fit to cos 9k of the form [3| 



-Fl cos^ Ok + t(1 - Fl){1 - cos^ 0k) 



(1) 



determines F^, the K* longitudinal polarization fraction. 
We also describe measurements of the distribution of the 
angle 9g between the £~^{i~) and the B{B) direction in 
the £^£^ rest frame. A fit to cos 9i of the form [3|] 

^Fl(1-cos2 0,) + |(1-Fl)(1+cos2 9i)+Afb cos9, (2) 
4 o 

determines Afb, the lepton forward-backward asymme- 
try. These measurements arc done in a low q^ region 
0.1 < q^ < 6.25 GcV^/c'*, and in a high q^ region above 
10.24 GeV^/c^. We remove the J/ip and i}){2S) resonances 
by vetoing events in the regions q^ = 6.25-10.24 GeV^/c^ 
and q^ = 12.96-14.06 GeV^/c" respectively. 

The SM predicts a distinctive variation of Afb arising 
from the interference between the different amplitudes. 
The expected SM dependence of ApB and Fl on q^ along 
with variations due to opposite-sign Wilson coefficients 
are shown in Fig. [3l At low g^, where C7" dominates, 
Afb is expected to be small with a zero-crossing point 
at (7^ ~ 4 GeV^/c^ [J, la, Ifil . There is an experimental con- 
straint on the magnitude of Cf^ coming from the branch- 
ing fraction for 6 ^ S7 [6|, |7[ , which corresponds to the 
limit q^ -^ 0. However, a reversal of the sign of C7" is 
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q q q q 

FIG. 1: Lowest-order Feynman diagrams for b — > st^l~ 



allowed. At high q^ , the product of Cg" and Cfg is ex- 
pected to give a large positive asymmetry. Right-handed 
weak currents have an opposite-sign Cg'^Cfp which would 
give a negative Afb ^^ high q^ . Contributions from non- 
SM processes can change the magnitudes and relative 
signs of Cf\ Cg" and Cfn, and may introduce complex 
phases between them [3|, [Sj. An experimental determi- 
nation of Fl is required to obtain a model-independent 
Afb result, and thus avoid drawing possibly incorrect 
inferences about new physics from our observations. 

We reconstruct signal events in six separate ffavor- 
specific final states containing an e^e~ or fi^ ^J.^ pair, 
and a K*{892) candidate reconstructed as /v+tt^, A'+tt^ 
or A'gTr"*" (or their charge conjugates). To understand 
combinatorial backgrounds we also reconstruct samples 
containing the same hadronic final states and e^fi"^ pairs, 
where no signal is expected because of lepton flavor con- 
servation. To understand backgrounds from hadrons (h) 
misidentified as muons, we similarly reconstruct samples 
containing h fjL^ pairs with no particle identiflcation re- 
quirement for the h^ . 

We use a dataset of 384 million BB pairs collected 
at the T(4S') resonance with the BaBAR detector [9| at 
the PEP-II asymmetric-energy e^e~ collider. Track- 
ing is provided by a flve-layer silicon vertex tracker 
and a 40-layer drift chamber in a 1.5 T magnetic held. 
We identify electrons with a CsI(Tl) electromagnetic 
calorimeter, muons with an instrumented magnetic flux 
return, and K^ using a detector of internally reflected 
Cherenkov light as well as ionization energy loss infor- 
mation. Charged tracks other than identified e, ^ and 
K candidates are treated as pions. Electrons (muons) 
are required to have momenta p > 0.3(0.7) GeV/c in the 
laboratory frame. We add photons to electrons when 
they are consistent with bremsstrahlung, and do not use 
electrons that arise from photon conversions to low-mass 
e+e~ pairs. Neutral A'^ — > tt+tt^ candidates are required 
to have an invariant mass consistent with the nominal K'-' 
mass [10|, and a flight distance from the e+e~ interac- 
tion point which is more than three times its uncertainty. 
Neutral pion candidates are formed from two photons 
with E^ > 50MeV, and an invariant mass between 115 
and 155MeV/c^. We require Ar*(892) candidates to have 
an invariant mass 0.82 < M{Kt:) < 0.97GeV/c^ 

B — > K*t^£~ decays are characterized by the kine- 
matic variables toes ~ ^s/4 — p*^ and A£' = E^ — 
Y^/2, where p*^ and Eg are the reconstructed B mo- 
mentum and energy in the center-of-mass (CM) frame, 
and ^/s is the total CM energy. We define a fit re- 
gion TOES > 5.2CeV/c^ with -0.07 < A£: < 0.04 
(-0.04 < AS < 0.04) GeV for e+e" [li+^i-) final 
states in the low q^ region, and —0.08 < /S.E < 0.05 
(-0.05 < AE < 0.05) GeV for high q^. We use the 
wider (narrower) AE windows to select the e fj,^ {h /j,^) 
background samples. 

The most significant background arises from random 



combinations of Icptons from semileptonic B and D de- 
cays. In BB events the leptons are kinematically corre- 
lated if they come from B -^ D^*^(.v, D -^ K^*^lv. Un- 
correlated backgrounds combine leptons from separate B 
decays or from continuum e+e~ -^ cc events. We sup- 
press these types of combinatorial background through 
the use of neural networks (NN) . For each final state we 
use four separate NN designed to suppress either con- 
tinuum or BB backgrounds in either the low or high q^ 
regions, and different selections of NN inputs are used 
depending on q^ bin (low, high), the identity of the lep- 
tons in the final state (e, /z), and the type of background 
(BB, continuum). Inputs include: 

• event thrust; 

• ratio of second-to-zeroth Fox- Wolfram mo- 
ments [11| : 

• mES and AE of the rest of the event (ROE), com- 
prising all charged tracks and neutral energy de- 
posits not used to reconstruct the signal candidate; 

• the magnitude of the total event transverse momen- 
tum, which is correlated with missing energy due to 
unreconstructed neutrinos in background semilep- 
tonic decays; 

• di-lepton system's distance of closest approach 
along the z-axis, and separately in the xy-plane, 
to the primary interaction point; 

• vertex probability of the signal candidate and, sep- 
arately, of the di-lepton system; 

• the cosines in the CM frame of the angle between 
the B candidate's momentum and the z axis, the 
angle between the event thrust axis and the z axis 
(^thrust), the angle between the ROE thrust axis 
and the z axis (^^^^^^gt)' ^^'^ the angle between 

^fh?u^st and ^thrust. 

There is also a background contribution in the signal 
region from B —f D{K*-k)tt decays, where both pions 
are misidentified. The misidentification rates for muons 
and electrons are ~ 2% and ~ 0.1%, respectively, so this 
background is only significant in the /^^/i~ final states. 
These events are vetoed if the invariant mass of the K*ti 
system is in the range 1.84-1.90 GeV/c^. 

We optimize the NN and AE selections for each final 
state in each q^ bin to give the best combined statisti- 
cal signal significance in the rriEs signal region toes > 
5.27GeV/c^ for the sum of all six final states. After all 
these selections have been applied, the final reconstruc- 
tion efficiencies and expected yields for signal events (cal- 
culated using world average branching fractions [7|), as 
well as expected yields for background events in the sig- 
nal region, are shown in Table HI 

For each q^ region, we combine events from all six final 
states and perform three successive unbinned maximum 
likelihood fits. Because of the relatively small number 
of signal candidates in each q^ region, a simultaneous fit 
over TTiEs , cos 9k and cos Oi is unlikely to converge and 



TABLE I: Signal efficiencies (%), and expected signal and 
background yields for tties > 5.27GeV/c'^, for low and high 
(f regions. 





Signa 


lEff. 


Signal Yield 


Bkgd. 


Yield 


Mode 


low 


high 


low high 


low 


high 


K+TT^fl+fl- 


1.6 


3.1 


1.0 1.8 


0.7 


3.8 


K°s7V+^,+ fl- 


3.6 


5.5 


3.0 4.5 


0.3 


1.4 


K+Tv-fi+^j.- 


4.5 


8.1 


5.5 9.6 


0.0 


3.1 


K+7r°e+e- 


4.6 


5.3 


2.8 3.1 


1.7 


2.4 


K°7r+e+e- 


7.0 


5.4 


5.9 4.4 


0.3 


1.4 


K+TT'e+e- 


8.6 


10.3 


10.5 12.2 


1.7 


2.4 


Total Yield 




28.6 35.8 


4.8 


14.5 



a sequential fitting procedure is required. We initially fit 
the TOES distribution using events with toes > 5.2 GeV/c^ 
to obtain the signal and background yields, Ns and 
Nb respectively. We use an ARGUS shape [ij] with a 
free shape parameter to describe the combinatorial back- 
ground in this fit. For the signal, we use a Gaussian shape 
with a mean toes = 5.2791 ± 0.0001 GeV/c^ and a ^ 
2.60 ± 0.03MeV/c^, which are determined from a fit to 
the vetoed charmonium samples. In this and subsequent 
fits we account for a small contribution from misidentified 
hadrons by subtracting the K*h^iif events, weighted by 
the probability for the h^ to be misidentified as a muon. 
We also account in all fits for charmonium events that 
escape the veto, and for mis-reconstructed signal events. 
We estimate contributions from non-resonant Ktt decays 
by fitting events outside the K* mass window in the range 
0.7 — 1.1 GeV/c^. We find no signal-like events that are 
not accounted for by the tails of the resonant mass distri- 
bution, and thus do not expect any significant contribu- 
tion from non-resonant events within the mass window. 

The second fit is to the cosine of the helicity an- 
gle of the K* decay, cosOk, for events with toes > 
5.27GeV/c'^. In this fit, the only free parameter is Fl, 
with the normalizations for signal and combinatorial 
background events taken from the initial toes fit- The 
background normalization is obtained by integrating, for 
'TIES > 5.27GeV/c^, the ARGUS shape resulting from 
the TO,ES fit. We model the cosOk shape of the combina- 
torial background using e^e^ and i-i^ i-J-^ events, as well 
as lepton-flavor violating e+^^ and yL^e~ events, in the 
5.20 < Toes < 5.27GeV/c^ sideband. The signal distri- 
bution given in equation (1) is folded with the detector 
acceptance as a function of cos Ok , which is obtained from 
simulated signal events. 

The final fit is to the cosine of the lepton helicity angle, 
cosOi, for events with toes > 5.27GeV/c^. The only free 
parameter in this fit is Afb , with the signal distribution 
given in equation (2) folded with the detector acceptance 
as a function of cosO(. In this fit, the value of F^ is fixed 
from the result of the second fit, and normalizations for 



TABLE II: Resuhs for the B -^ J/ipK* control samples. ABF 
are the differences between the measured branching fractions 
and the world average value [10|. The previously measured 
Fl = 0.56 ± 0.01 [III, and the expected Afb = 0. 



Mode 



ABF (10-=*) Fl 



Af 



K+Tv-fi+tJ,- 



Ls7 

K+n- 



+0.09 ±0.12 0.54 ±0.03 -0.04 ± 0.05 

±0.02 ±0.11 0.55 ±0.02 ±0.00 ± 0.05 

-0.03 ± 0.07 0.56 ± 0.02 -0.02 ± 0.02 

±0.16 ± 0.10 0.54 ± 0.03 ±0.02 ± 0.03 

±0.07 ±0.10 0.55 ±0.02 -0.02 ± 0.04 

±0.02 ± 0.07 0.56 ± 0.02 ±0.01 ± 0.02 



TABLE IV: Systematic errors on the measurements of Fl and 
Afb in the K*£^£~ samples. 



Source 


Fl 


A 


FB 


of Error 


low q^ 


high q^ 


low q 


high q^ 


rriBS fit yields 


0.001 


0.016 


0.003 


0.002 


Fl fit error 






0.025 


0.022 


Background shape 


0.011 


0.008 


0.017 


0.021 


Signal model 


0.036 


0.034 


0.030 


0.038 


Fit bias 


0.012 


0.020 


0.023 


0.052 


Mis-reconstructed signal 


0.010 


0.010 


0.020 


0.020 


Total 


0.041 


0.044 


0.052 


0.074 



TABLE III: Resuhs for the fits to the Ki+T and K't+T 
samples. Ns is the number of signal events in the tties fit. 
The quoted errors are statistical only. 



Decay 


e 


Ns Fl Afb 


Ki+r 


low 


26.0 ± 5.7 ±0.04t°;^| 




high 


26.5 ± 6.7 ±0.20lgi* 


K*i+r 


low 


27.2 ± 6.3 0.35 ±0.16 ±0.24t;j:^« 




high 


36.6 ± 9.6 0.7lto:22 +0-76J:E!:g 



signal and combinatorial background events are identical 
to those used in the second fit. We constrain the cosSf 
shape of the combinatorial background using the same 
sideband samples as for the cos 9k fit. The correlated 
leptons from B — > D'^*^£v, D — > K^*Hv give rise to an 
TOES-dependent peak in the combinatorial background at 
cos 9^ > 0.7, and we consider this correlation in our study 
of systematic errors. No such correlation is observed for 
cos 9k- 

We test our fits using the large sample of vetoed char- 
monium events. The branching fractions (BF) and K* 
polarization for B -^ J/ipK* are well known [10|, [l3| , and 
Afb is expected to be zero. The results of the fits to the 
six final states arc all consistent with expected values (see 
Table [ITl) ■ We further test our methodology by perform- 
ing the TOES and cos 9i fits on a sample of B^ -^ K'^£'^£~ 
decays. The results arc given in Table IIIII and are con- 
sistent with negligible forward-backward asymmetry, as 
expected in the SM and most new physics models [14| . 

We validate the fit model by performing ensembles of 
fits to datasets with events drawn from simulated signal 
and background event samples. The input SM values of 
Fl and Afb are reproduced with the expected statistical 
errors. A few percent of the fits do not converge due to 
small signal yields. We have also performed fits using 
signal events generated with widely varying values of Cf\ 
Cg" and CfJ covering the physically allowed regions of Fl 
and Afb, and find minimal bias in our fits. 

The systematic errors on the fitted values of Fl and 



Afb are summarized in Table IIVI The uncertainties 
in the fitted signal yields Ns, due to variations in the 
ARGUS shape in the toes fits, arc propagated into the 
angular fits. The errors on the fitted Fl values arc prop- 
agated into the Afb fits. We vary the combinatorial 
background shapes by dividing the sideband sample into 
two disjoint regions in tties- We vary the signal model 
usingsimulated events generated with different form fac- 
tors [a, |l5[, and with a range of values of Cj\ Cg" and 
C|q, to determine an average fit bias. Finally, the mod- 
eling of mis-reconstructed signal events is constrained by 
the fits to the charmonium samples (Table ITT)) , where it 
is the largest systematic uncertainty. 

The final fits to the K*£~^i~ samples are shown in 
Fig. [21 The results for Fl and Afb are given in Ta- 
ble IIIII and are shown in Fig. [3l In the low q^ region, 
where we expect Afb ^ —0.03 and Fl ^ 0.63 from 
the SM, we measure Afb = 0.24l°;23 ± 0.05 and Fl = 
0.35 ± 0.16 ± 0.04, where the first error is statistical and 
the second is systematic. In the high q^ region, the SM 
expectation is Afb ^ 0.38 and Fl ^ 0.40, and we mea- 
sure Afb = 0.76l^;g ± 0-07 and Fl = 0.7ltJ5:22 ± 0.04, 
with a signal yield of 36.6±9.6 events. Theoretical uncer- 
tainties on the expected SM Fl and Afb values are gen- 
erally difficult to characterize in the high q^ region, and 
although under better control for 1 < g^ < 6 GeV^/c^, the 
extension of our low q^ region below 1 GeV^/c* makes es- 
timates of uncertainties there difficult also. The quoted 
values are obtained using our implementation of the 
physics models described in [4|, [l5| , corresponding to the 
SM curves in Fig. H 

The expected SM value of Cj'J at next-to-next-to- 
Icading logarithmic (NNLL) order is Cl'^ = -4.43 [3]. 
A more recent NNLL calculation which evaluates con- 
tributions from the full set of seven form factors gives 
Qq = —4.13 [17|. The magnitude of possible contri- 
butions from new physics to Cio can be constrained if 
Afb > at high q^. By combining such a constraint on 
Afb with inclusive b — > s£~^£~ branching fraction results, 
an upper bound of |Cf^^| ^ 7 can be obtained, improving 
on an upper bound derived solely from branching frac- 




°\0S{B,f -1 



cos(e.) 



FIG. 2: K'e+r fits: (a) low q^ rn-ES, (b) high q^ tties, (c) 
low q^ cos 9k, (d) high g^ cos^k, (e) low q^ cos 9i, (f) high 
q^ cosOe; with combinatorial (dots) and peaking (long dash) 
background, signal (short dash) and total (solid) fit distribu- 
tions superimposed on the data points. 
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FIG. 3: Plots of our results for (a) Afb and (b) Fl for the 
decay B — > K*t^(.~ showing comparisons with SM (solid); 
Cf = -Cf (long dash); CfCf^ = -CfCfi (short dash); 
Cf = -Cf,CfCfo = -CfClo (dash-dot). Statistical and 
systematic errors are added in quadrature. Expected Fl val- 
ues integrated over each q^ region are also shown. The Fl 
curves with Cg"C"^^ = —Ctf^Cfo are nearly identical to the 
two curves shown. 



tion results of \C^^ 



, ^ 10 [iSl. O ur results are consistent 
with measurements by Belle 19[ , and replace the earlier 
BaBAR results in which only a lower limit on Afb was 
set in the low g^ region [20[. 
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